THE ABUNDANCE SPREAD IN THE BOOTES I DWARF 

SPHEROIDAL GALAXY 

John E. Norris^, Gerard Gilmore^, Rosemary F.G. Wyse^, Mark I. Wilkinson^'^, V. 
Belokurov^, N. Wyn Evans^, Daniel B. Zucker^ 

ABSTRACT 

We present medium-resolution spectra of 16 radial velocity red-giant mem- 
bers of the low-luminosity Bootes I dwarf spheroidal (dSph) galaxy that have 
sufficient S/N for abundance determination, based on the strength of the Ca II 
K line. Assuming [Ca/Fe] ~ 0.3, the abundance range in the sample is A[Fe/H] 
~ 1.7 dex, with one star having [Fe/H = -3.4. The dispersion is a"([Fe/H]) = 
0.45 ± 0.08 - similar to those of the Galaxy's more luminous dSph systems and 
uo Centauri. This suggests that the large mass (> 10'^ Mq) normally assumed 
to foster self-enrichment and the production of chemical abundance spreads was 
provided by the non-baryonic material in Bootes I. 



Subject headings: galaxies: dwarf — galaxies: individual (Bootes I) — galaxies: 
abundances 



1. INTRODUCTION 

The Galaxy's dSph satellites continue to play a challenging and pivotal role in our 
understanding of the ACDM paradigm and the manner in which the Milky Way formed. 
Not only do there appear to be too few of these systems in comparison with the predictions 
of CDM satellite dark halos (Klypin et al. 1999; Moore et al. 1999), but for those who see 
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the dSphs as possible building blocks for the Galaxy's halo (e.g. Carollo et al. 2007), their 
stellar populations pose problems. Not only are their age distributions too young (Unavane 
et al. 1996), their chemistry, at first blush, is wrong, at least as far as e.g. [a/Fe] is concerned 
(Venn et al. 2004, and references therein ). The report by Helmi et al. (2006) of a dearth of 
stars with [Fe/H] < -3.0 in the more luminous dSphs, if correct, offers no comfort to those 
who seek the origin of the now well-established Galactic extremely metal-poor stars having 
-4.0 <[Fe/H] < -3.0. 

The recent discovery of very low luminosity dSph galaxies, down to My ~ -4 (e.g. 
Belokurov et al. 2006, 2007), but with associated non-baryonic masses of 10^ Mq within the 
extent of the stellar distribution (e.g. Gilmore et al. 2007) hints there is still much to be 
understood before we have a complete understanding of the process of the formation of our 
Galaxy. 

An important feature of the more luminous of the Milky Way's dSphs is that there 
exists a well-established range in metallicity within each system, with dispersions of order 
0.3-0.5 dex (e.g. Helmi et al. 2006)0. The purpose of this Letter is to present evidence that 
similar behavior pertains to one of the less luminous systems, Bootes I, with My = -5.8 
(Belokurov et al. 2006). We shall argue that Bootes I exhibits a range in metallicity of ~ 1.7 
dex, (T([Fe/H]) = 0.45, and stars as metal-poor as [Fe/H] = -3.4. 



2. SPECTROGRAPHIC OBSERVATIONS AND DATA REDUCTION 

Spectra of ~ 330 Bootes I candidate members were obtained with the Anglo- Australian 
Telescope/ AAOmega fiber-fed spectrograph (]http : / / www. aao . gov . au/lo cal/ ' www / aaomega / ) 
combination during 2006 May 23-29 and 2007 April 18-20. In the blue, the spectra have 
resolution R = 5000, and cover the wavelength range 3850-4540 A, while in the red the 
corresponding numbers are 14000 and 8340-8840 A, respectively. 

Candidates were selected from the SDSS DR4 data set, using the color magnitude dia- 
gram selection mask illustrated in Figure 2 of Belokurov et al (2006). Stars in the magnitude 
range 17 < g < 20.5 were selected for observation with AAOmega, with a deliberate effort 
to select candidates up to one degree from the galaxy centre, four half-light radii, given 
the galaxy's apparently distorted and elliptical image. Only limited observational material 
was obtained in 2006. Deeper data, using essentially the same target selection, were ob- 



^Kirby et al. (2008) have very recently reported abundance spreads in several of the Galaxy's less luminous 
dSph systems. 
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tained in 2007. The 2006 run was the first major visitor use of the new AAOmega facility. 
These data sets were used to optimise and enhance the data reduction system, 2dfdr (see 
http://www.aao.gov.aU/AAO/2df/software.html#2dfdr ) , in collaboration with local exper- 
tise, leading to its current excellent performance. Final data calibration and reduction used 
what is now the public 2dfdr system. 

3. ANALYSIS 

3.1. Radial Velocities 

Heliocentric radial velocities were determined using the HCROSS routine of the FIGARO 
package (see http://www.starlink.rl.ac.uk/star/docs/sun86.htx/ node425.html) . This per- 
forms a cross-correlation between the program stellar spectrum and a template, to obtain 
the relative radial velocity. An associated confidence level and formal error are estimated 
and we accept only those with confidence = 1 (see Heavens 1993). We excised the strong 
Ca II H & K lines, although this made an insignificant difference in most cases. Two tem- 
plate spectra were used, of twilight sky and a G-giant star, both obtained during the same 
observing run as the Bootes I candidates. The velocities from the two templates were in 
general consistent within the errors and here we report those using twilight sky as template. 
We thus obtained velocities for 98 objects, for which our mean formal error is 3.6 km s~^. 
The internal accuracy on one observation, from repeat observations of 4 stars from 2006 and 
2007, is 7 km s~^ with a mean offset of —1.5 km s^^. Our external errors may be estimated 
by comparison with Martin et al. (2007). For 5 stars in common, observed by us in 2007, 
the mean offset is 4.6 km s~^, with a = 1.8 km s^^. 

3.2. Selection of Bootes I Members 

For the purposes of the present work we choose to consider objects for abundance 
analysis only if their spectra have net counts greater than ~ 200 per 0.34 A pixel at 4150 A 
and lie within the radial velocity range 90-115 km s"^. The first criterion is set to obtain 
a reliable abundance estimate based on the Ca II K 3933 A line, while the second admits 
objects within ± 1.8a of the galaxy's systemic velocity (we adopt a velocity dispersion of 
7 km s~^). Data for the 16 objects within our sample that satisfy these criteria are presented 
in Table 1, where columns (l)-(5) contain identification, coordinates, radial distance from 
galactic center, and heliocentric radial velocity, (6)-(7) present SDSS colors go and (g-r)o 
(we adopt E(B-V) = 0.02, following Belokurov et al. 2006), and column (8) contains the 
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net counts obtained per 0.34 A pixel at 4150 A. For four objects, data are available from 
both the 2006 and 2007 sessions - in these cases the second entry in the table refers to data 
obtained in 2006. 

We believe our sample has very small field contamination. In velocity windows at 
lower (65-90 km s~^) and higher (115-140 km s~^) values, and for our count restriction, we 
find a total of 4 objects (one clearly too strong-lined to be a member). Standard Galactic 
components of thin disk, thick disk and halo have predicted mean velocitie^ in this line-of- 
sight, for a typical distance of dwarf stars of 3 kpc, of -5, -25, and -130 km s~^, respectively, 
well-separated from Bootes I. 

3.3. Chemical Abundances from the Ca II K line 

For the putative Bootes I members in Table 1 we have measured the Ca II K line- 
strength index, K', and the G band index, G', defined by Beers et al. (1999), which we 
present in columns (9) and (10) of Table 1. Beers et al. provide a formalism, based on 
Galactic field stars and globular clusters that permits one to estimate [Fe/H] given K' and 
(B-V)o, well-calibrated down to [Fe/H] = -4.0, which we adopt here. The reader should 
be aware that the implicit assumption, that the same [Ca/Fe] vs [Fe/H] relationship applies 
within both the metal-poor Galactic objects used for the calibration and the dSph satellites, 
is not consistent with high-resolution abundance studies (e.g. Venn et al. 2004) at, say, 
the 0.3 dex level for the dSph systems studied to date. Here we shall cite [Fe/H] values 
obtained from the Beers et al. calibration, assuming [Ca/Fe] ~ 0.3; allowance for a putative 
spread in [Ca/Fe] in Bootes I would increase any inferred dispersion of [Fe/H]. To use the 
Beers et al. formulation we have obtained B-V values using the transformation (B-V)o = 
1.197x (g-r)o-|-0.049 valid on the range 0.70 < B-V < 1.3, determined from observations of 
the red giants in the Galactic globular clusters M13 and NGC 24lfl 

The resulting [Fe/H] values are presented in column (11) of Table 1. We find ([Fe/H]) = - 
2.51 ± 0.13, which agrees well with Munoz et al. (2006, [Fe/H] = -2.5), Martin et al. (2007, 



^See erratum in the appendix. 

•^Based on data kindly made available to us by H. L. Morrison (B-V and g-r photometry 
of Stetson (2008, http://wwwl.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/community/STETSON/standards/) and 
An et al. (2008), and cluster membership for NGC 2419 from Peterson (1985), Shetrone et 
al. (2001), and Suntzeff et al. (1988)). We prefer this calibration to those of Lupton (2005), 
http://www.sdss.org/dr5/algorithms/sdssUBVRITransform.html) and Zhao & Newberg (2006), from stars 
in the general field and metal-poor objects having (B-V) < 0.70, respectively. 



- 5 - 



[Fe/H] = -2.1), and Feltzing et al. (http: //www.mpa-garching.mpg.de~garcon08/program.html 



[Fe/H] = -2.7). Small sample statistics (Keeping 1962) of the abundance differences for the 
four stars in Table 1 observed in both 2006 and 2007 yields a standard deviation of a single 
observation of 0.19 dex0, while comparison of results for 7 stars in common with Feltzing et 
al. (assuming equal errors in their work and ours) yields 0.35 dex. Further consideration 
of Table 1 reveals the principal result of our investigation: among the 16 putative galaxy 
members there is a spread of A [Fe/H] ~ 1.7 - considerably larger than expected from the 
above measurement accuracies. The abundance dispersion for our 2007 data, corrected for 
measurement error of 0.19 (0.35) dex is 0.49 (0.40): we adopt (T[Fe/H] = 0.45 ± 0.08 for 
Bootes I. 

The large range in abundance can also be appreciated in Figure 1, which presents the 
spectra for 12 stars in Table 1. Figure 2 compares the most metal-poor object in our sample 
(Boo-1137, with [Fe/H] = -3.4) with classical metal-poor stars in the range -4.2 < [Fe/H] 
< -2.5. The low metallicity of Boo-1137 is well-confirmed|f] 



3.4. The Calcium Triplet 

To complement the Ca 11 K line data we present equivalent widths of the Ca 11 triplet 
8542 A line in column (12) of Table 1. (Poor sky subtraction prevents measurement of the 
triplet 8662 A line.) There is a clear and strong correlation between the K line and triplet 
line data in Table 1. The question that needs to be addressed is whether the triplet data 
support the K line abundances. 

To our knowledge, all abundance calibrations based on the triplet follow Armandroff 
& Da Costa (1991, hereafter AD) and are anchored to the abundances of Galactic globular 
clusters, none of which has [Fe/H] < -2.5. (This method also assumes [Ca/Fe] ~ 0.3.) 
Extrapolation is thus necessary when using the triplet for the ~ half of the objects in 



A second estimate of our internal accuracy comes from ^ 10,000 dwarfs observed during our 2006 session 
(and analyzed with similar techniques), which serendipitously included 16 objects with [Fe/H] < -1.5, and 
g > 18.0, from SSDS DR6 (see Lee et al. 2008). The dispersion of abundance differences between the two 
sub-samples (over 8 independent pointings) is 0.23 dex. 

^Boo-1137 lies 24' from the center of the system, well outside the half-light radius of 13' (but along the 
major axis of the galaxy's elongation). Given the rareness of objects with [Fe/H] < -3.0, however, it is very 
likely to be a member. From the HK survey for metal-poor stars. Beers et al.(1992) find 5 stars with [Fe/H] 
< -3.0, V < 14, and B-^V > 0.7, over 2300 deg^ For the HES, Christlieb (2008, priv. comm.) finds 9 stars 
with B < 14.5 and B-V > 0.7, over 7000 deg^. We estimate one should expect only ~ 0.02 halo giants 
having [Fe/H] < -3.0 within the central 30' of Bootes I. 
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Table 1 having [Fe/H] < -2.5. To make progress we proceed as follows. Based on spectra 
of standard globular cluster giants with [Fe/H] < -1.4 described by Norris et al. (1996), we 
transform W8542 to the AD triplet index Wad = W8542+W8662 using the relationship Wad 
= 2.34xW8542- The pivotal diagram of AD is their plot of Wad vs V-Vhb (V is the star's 
magnitude and Vhb refers to the horizontal branch of the cluster under consideration). To 
obtain V we transform the go and (g-r)o values in Table 1 following Lupton (2005, op.cit.) 
For Bootes I, Vhb = 19.55 (Dall'Ora et al. 2006). All but two of the objects in Table 1 
fall below the M15 ([Fe/H] = -2.2) data of AD (their Figure 3) in their uncalibrated region, 
reasonably consistent with only four stars in our Table 1 having [Fe/H] > -2.2. 

AD also define the reduced index W' = W8542+W8662 +0.619 x (V- Vhb), and determine 
the calibration [Fe/H] = [Fe/H]AD = 0.326W'-2.706, valid for 1.6 < W' < 4.2 (-2.2 < [Fe/H] 
< -1.3). Using these relationships for the data in Table 1, we obtain ([Fe/H])AD = -2.32 
and cr[Fe/H]AD = 0.20, somewhat different from the values obtained from the Ca II K line. 
For Boo-1137, with [Fe/H]ca 11 k = -3.4, we find [Fe/H]ca 11 triplet = -2.7. 

We found no model-atmosphere calibrations of the triplet in the literature for giants with 
[Fe/H] < -2.5. This is perhaps not surprising: in this regime the triplet is relatively weak, 
core-dominated, and challenging to model atmosphere analysis. Using the spectrum synthesis 



code MOOG of Sneden (2008, http: / /verdi.as.utexas.edu/moog.html[ ) , models of Kurucz 



(1993), atomic data from VALD (http://www.astro.uu.se/~vald/), and a range of van der 



Waals damping values for the triplet, together with T^s, log g, and colors from the Yale- 



Yonsei isochrones (http: / /www. astro. yale.edu/demarque/yyiso. html ) we computed spectra 



for [Fe/H] = -1.5 to -4.0, and plotted isoabundant loci in the Wad vs (V-Vhb) plane. The 
agreement between observation and theory leaves something to be desired: while the observed 
loci are essentially parallel, as are the theoretical ones, the latter have a shallower slope. Most 
importantly, the theoretical lines move closer together at lowest abundance. We use these 
results to indicate, approximately, the improvement we believe necessary to extrapolate the 
AD relationship. For [Fe/H] > -2.5, we adopt [Fe/H] = [Fe/H]AD = 0.326xW -2.706, and, 
for -4.0 < [Fe/H] < -2.5, [Fe/H] = [Fe/H]AD -1.03 + 0.83xW' -0.15xW'2. Then for the 
data in Table 1, we find ([Fe/H]) = -2.64, a[Fe/H] = 0.45, and [Fe/H] Boo-1137 = -3.8. 



4. DISCUSSION 

Despite its small baryonic mass, Bootes I has an abundance dispersion very comparable 
with those of the more massive dSphs (cr[Fe/H] = 0.3-0.6, e.g. Helmi et al. 2006) and 
the most massive (3x10^ M0) Galactic globular cluster u Cen ((T[Fe/H] = 0.3, Norris et 
al. 1996). With a current stellar mass of ~ 4x10^ Mq (Fellhauer et al. 2008), however. 
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the chemical inhomegeneity of Bootes I stands in stark contrast against the heavy-element 
homogeneity of Galactic globular clusters of similar stellar mass. The obvious solution is 
that Bootes I sits within the deep potential well of a dark halo mass - of order 10^^^ Mq. 
Given that chemical inhomogeneity within globular clusters (which contain no dark mass) 
appears to occur only in massive objects such as u Gen, it seems reasonable to assume that 
a requirement for the creation of heavy-element chemical inhomogeneity is that Bootes I 
should have a potential well deep enough to retain the ejecta of Type II supernovae for 
periods sufficiently long for them to be incorporated into subsequent stellar generations. 

There is a test of this: sclf-cnrichmcnt in the lowest luminosity systems will be greatly 
affected by sampling noise. As demonstrated by the dynamical modelling of FcUhaucr et 
al. (2008), if Bootes I ever had a dark matter halo, it still has it and has lost very httle of 
its stellar mass. That is to say, Bootes I has always had relatively little stellar material. If 
one normalizes abundance dispersion with respect to luminosity or current baryonic mass, 
the specific abundance dispersion of Bootes I is some 10-100 times larger than that of the 
more luminous dSph. Not many supernovae are needed to chemically enrich, say, 10^ M0 
of pristine material to an abundance of [Fe/H] = -2.5. One should expect very specific and 
different element patterns in the stars of Bootes I directly from this very small number of 
progenitors. One may be looking at the yields of individual supernovae in the abundances 
of these stars. 

The authors gratefully acknowledge the contributions of the AAOmega project team, 
in particular Rob Sharp, during this investigation. We thank H.L. Morrison for generous 
instruction in (B-V, g-r) - transformations. 

REFERENCES 

An, D. et al. 2008, astro-ph/0808.0001 (SDSS Paper 869) 

Armandroff, T.E. & Da Costa, G.S. 1991, AJ, 101, 1329 

Beers, T.C., Preston, G.W., & Shectman, S.A. 1992, AJ, 103, 1987 

Beers, T.C. et al. 1999, AJ, 117, 981 

Belokurov, V. et al. 2006, ApJ, 647, Llll 

Belokurov, V. et al. 2007, ApJ, 654, 897 

Carollo, D. et al. 2007, Nature, 450, 1020 

Cayrel, R. et al. 2004, A&A, 416, 1117 

DaU'Ora, M. et al. 2006, ApJ, 653, L109 



- 8 - 



Fellhauer, M. et al. 2008, MNRAS, 385, 1095 
Gilmore, G. et al. 2007, ApJ, 663, 948 
Heavens, A. 1993, MNRAS, 263, 735 
Helmi, A. et al. 2006, ApJ, 651, L121 

Keeping, E.S. 1962, Introduction to Statistical Inference (Princeton: Van Nostrand), 202 
Kirby, E.N. et al. 2008, ApJ, 685, L43 
Klypin, A. et al. 1999, ApJ, 522, 82 

Kurucz, R.L. 1993, CD-ROM 13, ATLAS9 Stellar Atmospheres Programs and 2 km/s Grid 
(Cambridge: SAO) 

Lee, Y.S. et al. 2008, astro-ph/0710.5645 

Martin, N.F. et al. 2007, MNRAS, 380, 281 

Moore, B. et al. 1999, ApJ, 522, L19 

Muiioz, R.R. et al. 2006, 650, L51 

Norris, J.E., Freeman, K.C., & Mighell. K.J. 1997, ApJ, 462, 241 
Peterson, R.C. 1985, ApJ, 297, 309 

Shetrone, M.D., Cote, P., & Sargent, W.L.W. 2001, ApJ, 548, 592 
Suntzeff, N.B., Kraft, R.P., & Kinman, T.D. 1988, AJ, 95, 91 
Unavane, M., Wyse, R.F.G., & Gilmore, G. 1996, MNRAS, 278, 727 
Venn, K.A. et al. 2004, AJ, 128, 1177 



Zhao, C. & Newberg, H.J. 2006, astro-ph/0612034 



This preprint was prepared with the AAS I^TJrjX macros v5.2. 



-9- 



A. Erratum 



An unfortunate confusion over coordinates resulted in incorrect predictions for the ve- 
locity distributions of foreground Galactic stars in §3.2. Instead of the quoted negative 
velocities, the correct predicted velocity distribution peaks close to zero heliocentric velocity 
for all Galactic populations, with only a small contribution (around 3% of the sample) in 
the velocity range of the Bootes I system. None of the analysis or conclusions of the paper 
is affected by this error. 
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Tablc 1. OBSERVATIONAL DATA AND DERIVED PARAMETERS FOR BOOTES I 



Star RA Dec r Vr go {g-r)o C4150 K' G' [Fe/H] W8552 

(2000) (2000) (') (kms-i) (A) (A) (A) 



(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


7 


13 


59 


35.53 


+14 20 23.7 


12.1 


104 


18.31 


0.70 


609 


7.49 


2.01 


-2.32 


1.35 


8 


13 


59 


38.62 


+14 19 15.9 


12.6 


109 


19.03 


0.56 


345 


4.53 


3.25 


-2.75 


0.84 


9 


13 


59 


48.81 


+14 19 42.9 


11.1 


108 


17.92 


0.80 


772 


6.62 


1.99 


-2.67 


1.16 


33 


14 


00 


11.73 


+14 25 01.4 


5.2 


105 


18.23 


0.74 


321 


5.17 


5.28 


-2.96 


1.20 


34 


14 


00 


21.11 


+14 17 27.9 


13.1 


104 


18.74 


0.65 


313 


4.94 


2.66 


-3.10 


1.20 


41 


14 


00 


25.83 


+14 26 07.6 


6.2 


98 


18.38 


0.70 


239 


8.13 


4.29 


-2.03 


2.17 


78 


14 


00 


14.73 


+14 13 13.9 


16.9 


102 


19.30 


0.64 


237 


6.10 


2.24 


-2.46 


1.14 




















202 


4.88 


1.48 


-3.02 


1.25 


94 


14 


00 


31.51 


+14 34 03.6 


7.4 


93 


17.53 


0.90 


538 


6.63 


0.14 


-2.79 


1.24 


117 


14 


00 


10.49 


+14 31 45.5 


2.1 


94 


18.19 


0.74 


283 


8.78 


3.46 


-1.72 


1.46 


121 


14 


00 


36.52 


+14 39 27.3 


12.0 


113 


17.92 


0.83 


237 


7.62 


3.23 


-2.37 


1.51 




















387 


7.43 


3.14 


-2.46 


1.54 


127 


14 


00 


14.57 


+14 35 52.7 


6.2 


99 


18.15 


0.76 


236 


9.47 


4.72 


-1.49 


1.67 




















346 


9.36 


3.96 


-1.49 


1.86 


130 


13 


59 


48.98 


+14 30 06.2 


4.1 


110 


18.21 


0.70 


337 


6.60 


2.13 


-2.55 


1.22 


911 


14 


00 


01.08 


+14 36 51.5 


7.0 


97 


17.94 


0.75 


237 


8.47 


3.44 


-1.98 


1.51 


980 


13 


59 


12.68 


+13 42 55.8 


48.8 


104 


18.51 


0.61 


396 


4.14 


0.82 


-3.09 


0.86 


1069 


13 


58 


53.22 


+14 06 57.8 


29.0 


111 


19.05 


0.56 


490 


5.57 


1.27 


-2.51 


0.96 




















200 


3.80 


1.52 


-2.91 




1137 


13 


58 


33.82 


+14 21 08.5 


24.0 


112 


18.11 


0.73 


859 


3.06 


1.09 


-3.45 


0.66 
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Fig. 1. — Spectra of 12 Bootes I red giants from Table 1 (continuum-normalized and broad- 
ened to resolution of FWHM = 2.5 A). Panels contain star name and go/(g-r)o/[Fe/H]. 
Objects with [Fe/H] < -3.0 are presented in blue, while red is used for those with [Fe/H] > 
-2.0. 
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Fig. 2. — Spectra of Boo-1137 ([Fe/H] = -3.4) and extremely metal-poor field giants of 
similar temperature (from ANU's 2.3m telescope; [Fe/H] and (B-V)o from Cayrel et al. 
(2004) and Beers et al. (1999), respectively). All comparison stars are hotter than Boo- 
1137, and would have slightly stronger Ca II K 3933 A lines than seen here, if they were as 
cool as it. 



